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Abstract 

Using a generalized procedure for obtaining the dispersion relation and the equation of motion 
for a propagating fermionic particle, we examine previous claims for a lightlike preferred axis 
embedded in the framework of Lorentz-invariance violation with preserved algebra. We show that, 
in a relatively high energy scale, the corresponding equation of motion is reduced to a conserving 
lepton number chiral equation previously predicted in the literature. Otherwise, in a relatively low 
energy scale, the equation is reduced to the usual Dirac equation for a free propagating fermionic 
particle. The new dynamics introduces some novel ingredients to the phenomenological analysis of 
the tritium beta decay. In particular, a modified cross section expression and the correspondent 
phenomenology of the end-point are evaluated. 
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Theoretical issues pointing to the possibility that many empirical successes of special 
relativity do need not demand Lorentz invariance of the underlying theory have been widely 
discussed in several frameworks [1-7]. The characteristic scale of such theories is likely to be 
associated with the Planck's energy E Pl ~ 10 19 GeV. However, the current attainable energies 
are minuscule compared to this scale, so that experimental signals are expected to be heavily 
suppressed, and no decisive evidence contradicting the exact Lorentz invariance has yet been 
experimentally detected. In this scenario, Cohen and Glashow pursue a different approach to 
the possible failure of Lorentz symmetry denominated very special relativity (VSR) [10, 11]. 
It is based on the hypothesis that the space-time symmetry group of nature is smaller than 
the Poincare group, consisting of space-time translations and one of certain subgroups of 
the Lorentz group. The formalism of VSR has been expanded for studying some peculiar 
aspects of neutrino physics with the VSR subgroup chosen to be the 4-parameter group 
SIM(2) [11]. Since neutrinos are now known to be massive, several mechanisms have been 
contrived to remedy the absence of neutrino mass in the Standard Model Lagrangian [12] . 
The framework of VSR admits the unconventional possibility of neutrino masses that neither 
violate lepton number nor require additional sterile states [11]. In the Dirac picture, lepton 
number is conserved with neutrinos acquiring mass via Yukawa couplings to sterile SU(2)- 
singlet neutrinos [13, 14]. In the Majorana picture, lepton number is violated and neutrino 
masses result from a seesaw mechanism involving heavy sterile states or via dimension-6 
operators resulting from unspecified new interactions [15, 16]. In spite of not being Lorentz 
invariant, the lepton number conserving neutrino masses are VSR invariant. There is no 
guarantee that neutrino masses have a VSR origin, but if so their sizes may be an indication 
of the magnitude of Lorentz-violating effects in other sectors, for instance, clS db SU. ggestion 
to the examination of the existence of a preferred axis in the cosmic radiation anisotropy 
[18]. 

In this Letter we report about an adequacy of the results of VSR into a Lorentz-invariance 
violation system reconstructed by means of modified conformal transformations acting on the 
Lorentz generators. By considering such frameworks allowing for the existence of a preferred 
frame, we derive expressions for the spectrum and the end-point of the beta decay, which 
can be used as an experimental probe of the peculiar way in which neutrinos experience 
Lorentz invariance. To reach this objective we combine each boost /rotation with an specific 
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transformation from which we introduce a preferential direction with the aid of a light-like 
vector defined as n M (= (1, 0, 0, 1)), n 2 = 0. The transformation has to be chosen as to bring 
an equation of motion which recovers the dynamics of the equation introduced in [11], 

(7%-^^)(l-7>(^) = 0, (1) 
V 2 p\n A J 

which admits the unconventional possibility of neutrino masses that neither violate lepton 
number nor require additional sterile states. We search for convenient unitary transforma- 
tions U acting on the usual Lorentz generators in order to recover the equation of motion 
for a free propagating fermionic particle. We expect that, in a relatively high energy scale, 
the corresponding equation will be reduced to the Glashow Eq.(l), and, in a relatively low 
energy scale, it will be reduced to the usual Dirac equation for a free propagating fermionic 
particle. 

Let us start with the definition of the momentum space M. as the four-dimensional vector 
space consisting of momentum vectors p^. The ordinary Lorentz generators act as 

L^x = p^d x - p\df, (2) 

where = d/dp 41 , and we assume the Minkowski metric signature and that all generators 
are anti-Hermitian (also /i, A = 0, 1, 2, 3, and i, j, k = 1, 2, 3 and the velocity of the 
light c = 1). The ordinary Lorentz algebra is constructed in terms of the usual rotations 
J 1 = e ijk L jk and boosts K i = L i0 as 

[J% K j ] = e ijk K k ; [J\ J 3 } = [K l ,K 3 ] = e l]k J k (3) 

In order to introduce the nonlinear action that modifies the ordinary Lorentz generators 
but, however, preserves the structure of the algebra, we suggest the following ansatz for a 
generalized transformation, 

D = (a(y) Plx + b(y)n^ (4) 

which acts on the momentum space as 

Dop^ = a{y)p^ + b{y)n^ (5) 

where y = p^n^. We assume the new action can be considered to be a nonstandard and 
nonlinear embedding of the Lorentz group into a modified conformal group which, as we 
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shall notice in the following for the case of main interest, despite the modifications, satisfies 
precisely the ordinary Lorentz algebra (3). To exponentiate the new action, we note that 

U = U(D) K* U(D) and f = U(D) J 1 U(D) (6) 

where the y-dependent transformation U (D) is given by U (D(y)) = exp[D(y)]. The nonlinear 
representation is then generated by U (D(y))op /Ji and, despite not being unitary (U (D(y))op^ ^ 
it has to preserve the structure of the algebra. Thus, when we assume 

[[L ltX ,D(y)],D(y)]=0 (7) 

we can reobtain a set of generators (in terms of k { and jj) which satisfy the ordinary Lorentz 
algebra of (3). At this point, in order to explicitly obtain the operator D(y) which satisfies 
the relation (7), we firstly compute the commuting relation 

[L mA , D(y)] = K ltX (a\y)p a + b'(y)n a )d a + b(y)d^, (8) 

for which we have defined the parameters k^x = p^rix — p\n^ and = ri\d^ —n^d\. From 
the above definitions we obtain the useful relations 

D k m a = a(y)Kn\, d^x D = a(y)d flX and D d^\ = 0, (9) 

which are essential in computing [[L^x, D{y)\, D{y)\. The first part of the r.h.s. of the Eq.(8) 
then leads to the commuting relation 

Kx(a'(y)px + b'(y)nx)dx, D{y)} = 

{ W 2 (y) - a(y) a"{y)) - a{y) a'(y)] p X d x 

+ [a'( y )(yb'(y)-b( y ))-ya(y)b"( y )]n X dx}, (10) 

and the second part gives 

[b(y)d^x, D(y)} = a(y)(b(y) -yb\y))d^- (11) 

In order to satisfy the condition for preserving the Lorentz algebra (7), the y-dependent 
coefficients can be obtained by evaluating the coupled ordinary differential equations: 

y{a'\y) - a(y) a"{y)) - a(y) a' (y) = (a.l) 

a\y){yb'{y) - b(y)) -ya{y)b"{y) = (a.2) 

a(y)(b(y)-yb'(y)) =0 (a.3) 

for which we have two types of solutions: 
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Type-I a(y) = and V b(y). 

Type-II a(y) = Ay n , n G 1Z and b(y) = By where A and B are constants with respective 
dimensions given by [[A]} = rrr n and [[£>]] = m" 1 . 

However, for Type-II solutions, it is difficult (and sometimes impossible) to constraint 
M (ay) to preserve the standard dispersion relation p 2 = m 2 . This objection makes the 
above solutions not so interesting as the Type-I solution. 

For a Type-I solution where a(y) = and b(y) = —am 2 /(I + lay), we can easily verify 
that 



a 



D(V) - =► D(y) o p, = _ T _^_ n, 

=> U(D(y)) o p M = p M(a) = Pjtt - — — — n M . (12) 

In spite of preserving the Lorentz algebra [17], these transformations clearly do not preserve 
the usual quadratic invariant in the momentum space. But there is a modified invariant 
\\U(D(y)) opuW 2 = M\a) which leads to the following dispersion relation, 

\\U(D(y)) o Pll \\ 2 =p 2 (a) =p 2 - = M\a) (13) 

Imposing the constraint p 2 = m 2 , which is also required by the VSR theory, we have the 
Casimir invariant 

M\a) = m2 . (14) 

for which the [/-invariance can be easily verified when we apply the transformation U [D(y)). 
The fact that the algebra of the symmetry group remains the same suggests that perhaps 
the standard spin connection formulation of relativity is still valid. In this sense, the above 
dispersion relation can also be obtained from the dynamic equation for a fermionic particle, 

[y (U(D(y)) op,) - M{a)\ u L (p v{a) ) = 

(771 " Q' \ 
7% " 1 + 2ay 7% " M («)J M^w) = 0. (15) 

Alternatively, as pointed out in [9] , for a comparative purpose to all these classes of Lorentz- 
violating models, dispersion relations may be derived from calculations in a theory such as 
loop quantum gravity [19]. 

By setting [[m]]" 1 values to a, for instance, a = ±l/m, ±m/e 2 pl (where e Pl is the Planck 
energy), we are able to analyze the low and the high energy limits. In the high energy limit 
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where ay » 1, the Eq. (15) is reduced to 



m 2 



7% - 2^7% " M («)J %(^(vsa)) = 0. (16) 

and since M\a) ~ ^t^t << """A i n spite of not being necessary, we can eliminate the 
dependence on a since the M(a) term becomes irrelevant as we can observe when we take 
the quadratic form of the above equation. Thus we recover the VSR Cohen-Glashow equation 
(1) [11] and its corresponding dispersion relation, as we have proposed from the initial part 
of this letter. In the low energy limit where ay « 1 (p u ( a ) — > Pu), the Eq. (15) is reduced 
to 

(l^Vii ~ m 2 ory^n^ - m) u L (p„) = 0, (17) 

whose the quadratic form is 

(p 2 - 2m 2 ay - m 2 ) u L (p„) w (p 2 - m 2 ) u L (p u ) = 0, (18) 

i.e. when ay « 1 the effective contribution from the second term of Eq. (17) is minimal 
and the equation can be reduced to the usual (low energy limit) Dirac equation for a free 
propagating particle. Such an important result could also be reproduced in a more direct 
way if we initially assumed a natural energy scale where ma << 1, for instance, when 
a = m/e 2 pl . Since M is reduced to m, the Eq. (17) is immediately reduced to the Dirac 
equation. 

Once we have established the novel dynamics, the measurement of /3-spectrum in the end- 
point region in tritium /3-decay is the classical method of direct determination of neutrino 
mass [21]. Since the dispersion relations is maintained, the differential of the decay rate for 
the d — > w e~ u e transition is related to the decay amplitude by [22] 

dr = g 2 £ \u(p e(a) W\i - 7r>(?v«))| 2 uXeawe,, 2 * 5 ^ ~ Ee ~ Ev) (19) 

spins e 

where E is the energy released to the lepton pair and G is the Fermi constant. For the novel 
dynamics, the phase space restriction is augmented by a change in the relevant matrix. We 
see that the weak leptonic charged current J M must be modified to ensure its conservation, 

a 2 m 2 , n^(n A 7 A ) 



J M = u(p e ) 



7 M + 



v(p v ) (20) 



2 (1 + 2ap*n x )(l + 2ap£n x ) 
The Lorentz violating a-dependent term leads to an entirely negligible effect near the end- 
point. It yields a maximal correction of order m u /m e (when am >> 1) which decreases as 
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the parameter a m diminishes. The first term, although not modified in form, differs from 
the standard current since the square matrix element v(p^)v(p^ a - ) ) is modified by 



v(P( a )MP( a )) 



-7s n 

2 Pis(a)7v 



(21) 



2 (l + 2aj£n A ). 

contrast to the standard expression. Now the summation over the spins in the Eq. (19) can 
be performed in order to give 



\u( Pe )j°(l-J 5 )v(P»)\ 2 = 8 



spins 



E. E, 



n a m L 



(l + 2ap>)nx) 



+ p e p l/ cos(6 e -9 v 



cos(0 e ) a ml 
(l + 2ap*n A ) 



(22) 

where 9 e ^ comes from the scalar product p e ^ ■ n = p e ^ cos(0 e ( I/ )). Substituting (22) into 
(19), the transition rates becomes 



dr 



16G 2 
(2tt) 6 



E. E-, 



n a m u 



+ Pe[pu cos(e e ~e u ) 



cos(6 e ) a ml 
(l + 2ap£n x )J_ 



x [(2vr(dcos(0 e ))p e2 dp e ) (27r(dcos(^))p 2 d^)] 6(E - E e - E v ) 



(23) 



where we have rewritten d 3 p e( v)/ (2E e ^) in terms of the spherical coordinates. By performing 
the angular integrations and the E u (p u ) integration, after some mathematical manipulations, 
we obtain 

~\ + 2u(E v +p v 



dT 

dp 



= Cp 2 {E uPu 



m: 



In 



2 [l + 2a{E u -p u ) 



(24) 



where we have suppressed the electron index e from the kinematical variables. In the above 
equation, E v = E -E = (K max +m e )-(K+m e ), p v = a/ (K max - K) 2 - m 2 and the constant 
C is defined as G 2 /2n 3 . The Kurie plot rate p~ 1 (dr/dp e ) as a function of the neutrino 
energy (E — E ) near the end-point of the tritium beta decay spectrum (K max = 18.6 keV) 
for m u = 1 eV is shown in the first plot of the Fig.l. Since the final state neutrinos are 
not detected in the tritium /3-decays experiments, for the electron spectrum, we get the 
incoherent sum 

^ = t\U ej \ 2 ^ 1 (25) 
dp ^-f 1 Jl dp y ' 

3 = 1 

for a superposition of two neutrino mass eigenstates with m 1 — 1 eV and m 1 — 2 eV[24], with 
mixing angle equal to 7r/4, in the second plot of the same figure. It is important to notice 
that, in a realistic phenomeno logical analysis, we can vary continuously the parameter a 
from to oo (VSR dynamics limit with am » 1), as well as add the necessary number 
of neutrino species in (25). In fact, the knowledge of neutrino mass spectrum is decisive 
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FIG. 1: The Kurie plot rate p~ 1 (dT/dp e ) as a function of the neutrino energy (E — E ) near the 
end-point of the tritium beta decay spectrum for m u = 1 eV in the first plot, and for m ul = 1 eV 
and m y2 = 2eV in the second plot. The dotted line is for a massless neutrino. The parameter a 
varies continuously from zero (the low energy or Dirac equation limit) to infinity (the high energy 
or the VSR limit). Here is defined as G/(2vr 3 ) 1 / 2 

for the understanding of the origin of neutrino masses and mixing. If in the KATRIN [23] 
experiment, which is under preparation at present, a positive effect due to the neutrino 
mass will be observed, we will have Tn u fp\ ~ m Ul23 [25]. The KATRIN experiment and its 
predecessors measured the integrate energy spectrum from the end-point downward which 
is proportional to 

fK max —m v j-p 

T(K) = j K lK dK < 26 > 

where K(E U = E — E ) is implicitly denned above as the the electron kinetic energy K = 

8 
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FIG. 2: The integrated end-point difference T mv=0 (K) — T(K) as a function of the neutrino energy 
{—E v = K — K max = E — Eq) near the end-point of the tritium beta decay spectrum for m u = 1 eV 
in the first plot, and for m ul = leV and m u2 = 2eV in the second plot. Again we vary the 
parameter a from zero to infinity in order to obtain the fine-tuning correction to integrated energy 
spectrum from the end-point downward. Here C is defined as G 2 /2tt s 

E — m e = -Kmax + E u . By observing that K dK = E dE = p dp, we can numerically evaluate 
the above integral in order to obtain the Fig. 2. The effect of neutrino mass is conveniently 
expressed as the difference from the massless case in terms of T m ^ =0 (K) — T(K) as a function 
of the neutrino energy (E — E ). By knowing the experimental inputs, we could fit the 
results for the best values of a and m and compare with the standard result (a = 0). 

To summarize, we have examined previous theoretical claims for a preferred axis at 
?2 M (= (1,0,0, 1)), n 2 = in the framework of Lorentz invariance violation by generalizing 
the procedure for obtaining the equation of motion for a propagating fermionic particle. We 



a = m = 1 eV 
a = 0.1 

a = 1 
a» 1 
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have shown that, in a relatively high energy scale, the corresponding equation of motion is 
reduced to a conserving lepton number chiral equation previously predicted in the literature 
[11], and, in a relatively low energy scale, it is reduced to the usual Dirac equation for a free 
propagating fermionic particle. For fermionic particles with the dynamics described by the 
Eq. (9), the phase space kinematic restriction is modified by a change in the relevant matrix 
element. In order to incorporate some discrete spatial and causal structures at the Planck 
energy scale, the action which leads to Lorentz invariance with an invariant energy scale (e Pl 
or m 2 /e Pl ) can be taken into account simultaneously with the action here proposed. Ideally, 
the formalism we are discussing, in the sense analogous to that of the VSR, can be used to 
compare experiment and theory, as well as to extrapolate between predictions of different 
experimental measurements. For this aim, we have quantified the eventual modifications to 
which could be observed in the next generation of tritium beta decay end-point experiments. 
From the fit of the experimental data, it should be interesting obtaining an explicit relation 
between our free parameter a and the Planck energy. Moreover, the formal construction 
presented in this letter can be used for parameterizing two other possible phenomenologically 
observable modifications to neutrino physics: (i) eventual modifications to the predictions 
to neutrinoless double beta decay [20], and (ii) small modifications to the oscillation picture 
due to Lorentz-violating interactions that couple with the active neutrinos and eventually 
allow the complete explanation of neutrino data. With the compelling evidence for massive 
neutrinos from recent oscillation experiments, one of the most fundamental tasks of particle 
physics over the next years will be the determination of the absolute mass scale of neutrinos. 
Beside of being part of a phenomenology of quantum gravity effects, as opposed to directly 
having a fundamental significance, this formalism may eventually contribute for determining 
the absolute value of neutrino masses, which will have crucial implications for cosmology, 
astrophysics and particle physics. 
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